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Testing the Nature of the D*j(2317) + and £> sJ (2463) + States Using Radiative Transitions 
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The Babar and CLEO collaborations have recently observed states decaying to Dfir and D* + 7r° 
respectively and suggest the possible explanation that they are the missing P-wave cs states with 
J = + and 1 + . In this note we compare the properties of the D*j(2317) + and Z? s j(2463) + states 
to those expected of the cs D* and D s \ states. We expect the D* and D s i with the reported 
masses to be extremely narrow, V ~ O(10 keV), with large branching ratios to D*7 for the D*o and 
to -D*7 and D s 7 for the D s \. Crucial to this interpretation of the Babar and CLEO observations 
is the measurement of the radiative transitions. We note that it may be possible to observe the 
D s i(2536) in radiative transitions to the D*. 

PACS numbers: 12.39.-x, 13.20.Fc, 13.25.Ft, 14.40.Lb 



I. INTRODUCTION 

Over the last decade there has been considerable 
progress in our understanding of mesons, strongly inter- 
acting bound states of quarks and antiquarks. Mesons 
made of one heavy and one light quark have played an 
important role However the theoretical predictions 
have not been sufficiently tested by experimental data to 
say that we truly understand the strong interaction. This 
situation has recently been highlighted by the discovery 
of astate, the D* J (2317)+, with mass 2317 MeV decaying 
to Dfir by the Babar Collaboration at the Stanford Lin- 
ear Accelerator Center (SLAC) Q and a second state, the 
L> sJ (2463)+, with mass 2.463 GeV decaying to L»*+7r° by 
the CLEO Collaboration at the Cornell Electron Storage 
Ring 0- These states have also been observed by the 
Belle Collaboration at KEK 0. The J D* J (2317)+ was ob- 
served in the inclusive Dfir invariant mass distribution 
,2j . The state has natural spin-parity and the Babar col- 
laboration suggest it to be a J p — + state based on its 
low mass. The quantum numbers of the final state indi- 
cate that the decay violates isospin conservation. Babar 
found no evidence for the decay D*j(2317) + — * D* s + "f or 
Dfjj and although they found no evidence for the de- 
cay D*j(2317)+ -> £>+7 they see a peak near 2.46 GeV 
in the D+ir°j mass distribution but do not claim this as 
evidence for a new state. 

The CLEO collaboration subsequently reported on a 
signal in the D* + tt° channel at a mass of 2.463 GeV 
which they refer to as the D s j(2463) + Because the 
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D s j(2463) lies above the kinematic threshold to decay to 
DK but not D*K the narrow width suggests the decay 
to DK does not occur. Since angular momentum and 
parity conservation forbids a 1 + state from decaying to 
two pseudoscalars, CLEO suggests the compatability of 
the £> sJ (2463) with the J p = 1+ hypothesis. CLEO puts 
limits on the widths of the L»* (2317)+ and D s , 7 (2463)+ 
of T < 7 MeV at 90% C.L. 0. More importantly for 
the purpose of this analysis, they give limits on radiative 
transitions of the L>* 7 (2317) + and £> sJ (2463)+ to D*+j 
and -D+7 final states. 

The simplest interpretation is to identify these states 
as the missing j — 1/2 members of the cs L = 1 mul- 
tiplet where j — 1/2 is the angular momentum of the 
s-quark. The observation of these states is surprising 
because they are narrower than expected, are observed 
in isospin violating D s ir° and D*ir° final states, and are 
lower in mass than expected by most (but not all) cal- 
culations IllEEIllilllllimHlIl- The Babar 
and CLEO observations have led to conflicting interpre- 
tations. Although the observed mass for the D s o candi- 
date (the J p = + member of the ground state L = 1 
cs multiplet) is consistent with some predictions of chi- 
ral quark models 0, m which broken chiral symme- 
try views them as the positive-parity partners of the D s 
and D* states, it is considerab ly lo wer than expected by 
most quark models H, IE IE IS El an d lattice QCD cal- 
culatio ns fl 2 | IT3 | . This has led to considerable interest 
[illll7l liall'9| including the proposal that the L»j(2317)+ 
is a multiquark state [IE , possibly a DK molecule 
analogous to the KK interpretation of the /o(980) and 
a (980). 

In this letter we confront the cs L = 1 D* and D s \ 
interpretations of these states with the theoretical expec- 
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tations for conventional cs states. 



II. SPECTROSCOPY 

Mass predictions are an important test of QCD mo- 
tivated potential models as well as other calculational 
app roaches for hadron spectroscopy 0, IE S IS S \Ei EH 
In QCD-motivated potential models the spin- 
dependent splittings test the Lorentz nature of the con- 
fining potential with different combinations of Lorentz 
scalar, vector,. . . interactions H H SH H IH El • Fur- 
thermore, the observation of heavy-light mesons is an im- 
portant validation of heavy quark effective theory |22ll23j | 
and lattice QCD calculations 0, El- In Table I we 
summarize predictions for the P-wave cs states. The 
two J — 1 states are linear combinations of 3 Pi and 
1 Pi because for unequal mass quarks, C is no longer a 
good quantum number. 
D\. Most, but not all, 
the D* and the missing £q state to be substantially 



TABLE I: Predictions for the P-wave cs states. The J = 1 
states are linear combinations of the 3 Pi and 1 Pi states. In 
column 3 we list D\ which we take to be the higher mass state 
of the two J = 1 physical states and D[ the lower. PDG refers 
to the particle data group |24| and LGT refers to the lattice 
gauge theory result. 



Reference 



3 P D\ D{ 



We label these as the D\ and 
models predict the masses of 



Babar [2] 
CLEO [3] 
PDG [24] 
GI [5, 6] 
ZVR [7] 
EGF [8] 
DE [9] 
GJ [10] 
LNR fit] 

LW [LGT] [12] 2.499 2.511 2.500 2.554 
GB [LGT] [13] 2.437 



2.32 

2.463 
2.535 2.574 

2.48 2.56 2.55 2.59 

2.38 2.52 2.51 2.58 

2.508 2.569 2.515 2.560 

2.487 2.605 2.535 2.581 

2.388 2.536 2.521 2.573 

2.455 2.522 2.502 2.586 



higher than the masses reported by Babar |^ and CLEO gtate The fiabar collaboration observe d the L>%(2317) 
|3j| ■ Although it is possible that these models need re- 
vision it seems unlikely that they would disagree with 
experiment to such a large degree given their general 
success in describing the meson spectrum. A more se- 
rious problem is the large discrepancy with lattice QCD 
calculations which give M(D* n ) = 2499(13) (5) MeV and 

m 



If the D* sJ (2317) + and 



M(D s i)=2500(16)(2) MeV 
D s j(2463) + are identified as the missing a Po(cs) and 
Pi(cs) states it would pose a serious challenge for the 
lattice calculations. 

Quark model calculations [6| and heavy quark symme- 
try [22] predict that the 4 L = 1 cs mesons are grouped 
into two doublets with properties characterized by the 
angular momentum of the lightest quark, j = 1/2 and 
j = 3/2. The j = 3/2 states are identified with the pre- 
viously observed P 5 i(2536) ± and P sJ (2573) ± states 
while the j = 1/2 have not previously been observed. 
The j = 3/2 states are predicted to be relatively narrow 
0, in agreement with experiment [24|. In contrast, as- 
suming the higher masses predicted by the quark model, 
the j = 1/2 states are expected to be rather broad, de- 
caying to DK and D*K respectively with large ,5- wave 
widths. The large width is presumed to explain why 
they have yet to be observed. However, if these states 
are identified with the recently observed D*j(2317) and 
P s ,/(2463) states their masses would be below the DK 
and D*K thresholds so that they would be quite narrow, 
especially for mesons with such high mass. 



III. RADIATIVE TRANSITIONS 

While masses are one test of models of hadrons, tran- 
sitions probe the internal structure of the state. Com- 
parison between theory and experiment of the branch- 
ing ratios is an important test of any assignment for a 



in the D s ir° final state and report no observation of its 
decay via radiative transitions Q. The CLEO collabo- 
ration put limits on branching ratios of various radiative 
decays of the D*j(2317) + and P sJ (2463)+ relative to 
r(D; j (2317)+ -> D+n°) and r(P sJ (2463)+ -> D*+tt°) 
respectively 3]. Because the D* sJ (2317) + 's mass is below 
the kinematic threshold for the decay D s q — > DK, the 
only kinematically allowed strong decay is D s $ — > D s ir°. 
Likewise, the P s j(2463) is kinematically forbidden to de- 
cay to its expected dominant decay mode D s \ — > D* K 
so that the D s \ — > D*ir° is expected to be dominant . 
In both cases the decays D* Q — > D s tt° and D s \ — > D*tt° 
violate isospin and are expected to have quite small par- 
tial widths. Thus, the radiative transitions D* — > -D*7, 
P s i — > P*7 and D s i — > D s j would be expected to have 
prominent branching ratios. 

The PI radiative transitions are given by 



T(i -» / + 7 ) 
4 
27 



(1) 

-n , g ) 2 ^ (2J / + 1) \( 2s+1 Sj,\r\^Pj)\ 2 S lf 



where Sif is a statistical factor with S,, 



if 



1 for the tran- 



sitions between spin-triplet states, d[*J (IP) — > P*7 and 
D*(2S) — » D S ,/(1P)7, and = 3 for the transition be- 
tween spin-singlet states, D s \ — > P s 7, (cq) is an effective 
quark charge given by 



(eg) 



m s e c 



m c e- s 



(2) 



where e c = 2/3 and e s — 1/3 are the charges of the 
c-quark and s-antiquark given in units of \e\, m c = 
1.628 GeV, m s = 0.419 GeV are the mass of the c and 
s quarks taken from Ref. [5|, ot = 1/137.036 is the fine- 
structure constant, and lo is the photon's energy. The 
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matrix elements (S\r\P), given in Table II, were evalu- 
ated using the wavefunctions of Ref. Relativistic 
corrections are included in the El transition via Siegert's 
theorem pM l26l l27| by including spin dependent inter- 
actions in the Hamiltonian used to calculate the meson 
masses and wavefunctions. To calculate the appropriate 
photon energies the PDG [24| values were used for ob- 
served mesons while the predictions from Ref. [f| were 
used for unobserved states with the following modifica- 
tion. While splittings between cs states predicted by Ref. 
p| are in good agreement with experiment the masses are 
slightly higher than observed so to give a more reliable 
estimate of phase space, the masses used in Table II have 
been adjusted lower by 18 MeV from the predictions of 
Ref. Q. For the D* and D s \ states we give one set 
of predictions using the Babar and CLEO masses and 
a second set of predictions using the quark model mass 
predictions of Ref. p| . 

A final subtlety is that the J = 1 states are linear 
combinations of 3 Pi and 1 P\ because for unequal mass 
quarks, C is no longer a good quantum number. Thus, 

D 3 / 2 = l P l om6+ 1i Pywa.0 

D][ 2 = - 1 P 1 sin0+ 3 Picosl? (3) 

we use 9 = —38° and the conventions of Ref. [6j in 
calculating the widths in Table II which include the ap- 
propriate factors of cos 2 9 and sin 2 9 as appropriate. The 
resulting widths are given in Table II. 

In addition to the El transitions the Ml transitions 
D s \ — > D* "f can also take place. However, we found 
these partial widths to be quite small and unlikely to be 
observable. 



IV. STRONG TRANSITIONS 

The transition D* — > D s ir° is expected to be quite 
small as it violates isospin. Athough there are a num- 
ber of theoretical predictions for hadronic transitions be- 
tween quarkonium levels |2^, |3(3, I2J |U H^l we know 
of none for the transition D* — > D s ir°. To estimate this 
partial width we turn to known transitions and use exist- 
ing theoretical calculations for guidance. This approach 
should at least help us gauge the relative importance 
of this partial width. The only measured transition is 
ip{2S) -v J/ip(lS) + n with£ = 9.7xl0" 4 implying 
r(-0' -> J/-07T ) = 0.27 keV. A limit exists on the tran- 
sition T(2S*) T(15) + tt° of B(T(2S) -» T(1S)tt°) < 
1.1 x 10" 3 90% C.L. implying r(T(25) -> T(1S>°) < 
0.05 keV 0. The BR for the transition D* s -> D s +tt° is 
5.8 ± 2.5% but the total width is not known. We can es- 
timate the width by using the measured branching ratio 
B(D* -> D sl ) = (94.2±2.5)% with a quark model calcu- 
lation of the radiative transition DZ — > D s j. Combining 
the partial width given by Ref. @ of T(D* -> D sl ) = 
0.125 keV with the measured branching ratio 01 gives 
T(D* — > D s ir a ) ~ 7.7 eV. For comparison Goity and 



Roberts obtain T(D* S -> P s7 ) = 0.165 keV giving 
T(D* S -> D s tt°) = 10 eV (for the K = 0.45 solution) and 
Ebert et al. || find T(D* -> D s -y) = 0.19 keV giving 
T{D* D s n°) = 12 eV. 

For our first attempt to estimate T(_D* — > D s tt°) we 
rescale T(_D* — > D s tt°) assuming a fc 3 dependence for 
the partial widths and find r(L»; o -> D s tt°) ~ 2 keV. 
One should take this estimate with a grain of salt as 
the D* — > £> s 7r° is an 5 — > 5 transition with the final 
states in a relative P-wave while the D* Q — > Z? s 7r° tran- 
sition is a P — > 5 transition with the final states in a 
relative S'-wave so there are wavefunction effects we have 
totally ignored in addition to a generally cavalier atti- 
tude to kinematic factors. All we have attempted to do 
is establish the order of magnitude. 

A more relevant starting point is the transition 
/i c ( 1 Pi) — > J/-07T which is a P — > 5 spin-flip transi- 
tion which proceeds via the PI — Ml interference term 
in a multipole expansion of the gluonic fields, similar 
to the 3 Po —i- 1 Sq transition we are attempting to es- 
timate. Ko estimates T(h c —* J/tpir°) ~ 2.5 keV [33J. 
This transition is related to the transition tp' — » h c Tr° 
H3 HI for which Ko [H obtains B{^' -> /i c tt ) = 
3 x 10 -3 . For comparison Voloshin |2g finds B(ip' -» 
/i c 7r°) = 10~ 3 so that we should assume a factor of 3 
in uncertainty. These transitions are proportional to 
the pion momentum so that by rescaling the estimate 
of T(h c -> J/Vtt ) we find T{D* s0 -> P s vr°) ~ 2 keV. 
There are two important uncertainties in this estimate. 
The first is that the matrix elements are proportional 
to (S\r\P). Using the wavefunctions of Ref. |5j we 
find (l 3 P |r|l 1 5 ) cs /(l 3 5 1 |r|l 1 P 1 ) cc = 1.1. The sec- 
ond uncertainty is that the matrix elements are 0(a s ) 
so that the ratio of the widths go like (a s (cs)) / a s (cc)) 2 
which, given that the relevant energy scale is the light 
quark mass, could contribute an additional factor of 4 in 
the width. Given these uncertainties wc estimate that 
r(D* (2.32) -> P S 7T°) - 10 keV. We expect similar rates 
for the decays D s i — > D*tt° |3J|. In addition to the one- 
pion decay modes, the D s \ state can decay via two-pion 
transitions to the D s state. (The decay D* Q — > D s itit is 
forbidden by parity conservation.) Using Ko's estimate 
of the ratio T(h c -> J/tp + mr)/T{h c -> J/t/i + tt") ~ 0.16 
[33| we estimate T(D S \ — > D s mr) ~ 1.6 keV. The result- 
ing partial widths and branching ratios are summarized 
in Table II. 

For comparison we also include in Table II the partial 
widths and branching ratios expected for the l 3 p)(cs) 

1/2 

state with mass 2.466 MeV and the D sl state with mass 
2.536 MeV. The dominant decays for these masses are 
D* -> DK and D]{ 2 -> D*K with large partial widths. 
Although there is considerable uncertainty in the esti- 
mate of these widths 0, |3f| we do expect the P* and 

1 12 

D sl ~ states with these masses to be rather broad with 
small branching ratios for the radiative transition. These 
decays are S-wave so the widths scales linearly with the 
decay products momentum. 
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For completeness we also include in Table II other El 
transitions involving the cs P-wave states. We note that 
the D s i(2536) ± should have a relatively large branching 
ratio for its radiative transition to so that it may 

be possible to observe the D s i(2536) ± in this mode. 

CLEO [3j has obtained 90% C.L. limits on radiative 
transitions of the D*j(2317) and D s j(2463) which we 
summarize along with our predictions in Table III. 

V. OTHER POSSIBILITIES 

If the radiative transitions are not observed with BR's 
consistent with those of the conventional D* and D s i 
states what are the alternatives? One possibility sug- 
gested by the Babar collaboration is that the D* 7 (2317) + 
is some sort of multiquark state, either a DK molecule 
or a cqqs multiquark object. This seems to be a likely 
possibility which has much in common with the descrip- 
tion of the /o(980) and do (980) as multiquark states: The 
D* sJ (2317) + lies just below the DK threshold while the 
/o(980)/a (980) lie just below the KK threshold and 
both couple strongly to these nearby channels. This 
explanation has been promoted by Barnes, Close and 
Lipkin poj and is supported by a recent dynamical cal- 
culation by van Beveren and Rupp |2l) . Likewise, the 
D s j(2463) could be a D*K bound state similar to the 
K*K molecule interpretation advocated as the solution 
to the longstanding E/i puzzle |36| . 



The Babar and CLEO collaborations believe that they 
may be the missing J p — + and 1 + members of the 
L = l(cs) multiplet. However their masses are signifi- 
cantly lower than expected by most models and also lat- 
tice QCD calculations and would pose a serious challenge 
to these calculations. It is therefore important to test 
these assignments. If the D* 7 (2317) + and D sJ (2463) are 

conventional D* and D X J^ (cs) states we have argued 
that they should have very small total widths, O(10) keV, 
with large branching ratios to D*j (and D s j for the 

1 /2 

D sl ). It is therefore important to make a better deter- 
mination of the total width of these states and to search 
for the radiative transitions. In contrast, the absence of 
the radiative transitions and a relatively large total width 
of O(MeV) would support the D^K molecule designa- 

1 /2 

tions. In this case the conventional D* Q and D x ' states 
have yet to be discovered, presumably due to their large 
width. However, observation of their non-strange part- 
ners by the Belle collaboration [37J with their expected 
properties leads us to be hopeful that they can be found. 
Acknowledgments 



VI. CONCLUSIONS 

The discovery of the D* sJ (2317)+ and D sJ (2463) has 
presented an interesting puzzle to meson spectroscopists. 
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TABLE II: Predictions for partial widths and branching ratios for El transitions 2S — » IP and IP — > 2S 1 and strong decays in 

1/2 

the D s meson sector. For the D* and D sl states we show results for two sets of assumptions. In the first we associate the 
newly observed D* sJ (2.317) and D s j (2.463) with the D* and D^ 2 while in the second we show partial widths using the quark 
model predictions for these states's masses. For decays involving the D s i states we include the appropriate cos 2 8 and sin 2 6 
factors corresponding to eqn. 3 in the partial widths. The widths are given in keV unless otherwise noted. The masses come 
from the PDG |24j| unless otherwise noted. 
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"From Babar Ref.0. 

'Obtained using the wavefunctions generated from Ref. |jj 

c Masses taken from Ref. 5] with the modification that the pre- 
dictions have been adjusted downward by 18 MeV to give better 
agreement with the measured masses. The masses in Ref. |^ were 
rounded to 10 MeV. Here we round to 1 MeV. 

''Obtained by rescaling the result of Ref. 0] by phase space. 

Trom CLEO Ref.0. 

^Based on the PDG total width for the D sJ (2573) ± [H 
"The PDG gives Y < 2.3 MeV 90% C.L.. We used the width given 
by Ref. rescaled for phase space. 



TABLE III: Comparison of 90% C.L. limits on radiative tran- 
sitions obtained by CLEO with the predictions given in Ta- 
ble II. The BR's are with respect to the decay D* (2317) — > 
D s 7r° for the D*j(2317) and with respect to the decay 

D^ 2 (2463) -> D*w° for the P> sJ (2463). 



Transition 




Predicted 


CLEO [3] 


#1/(2317) 




0.19 


< 0.059 


D* SJ (2317) 


->D+>y 





< 0.052 


P> s j(2463) 


-» D*+ 7 


0.55 


< 0.16 


P> s j(2463) 




0.62 


< 0.49 


D s j(2463) 


-» D* sj (2317)7 1.2 x 10 -3 


< 0.58 
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